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ABSTRACT 

We report on numerical calculations of nonadiabatic eigenvalues and eigenfunctions for 
g-modes in ZZ Ceti variables. The spectrum of overstable t = 1 modes delineates the instability 
strip. Its blue edge occurs where ojt c « 1 for the n = 1 mode. Here u) is radian frequency and r c 
is about four times the thermal timescale at the bottom of the surface convection zone. As a ZZ 
Ceti cools, its convection zone deepens, longer period modes become overstable, but the critical 
CD , value of lot c separating overstable and damped modes rises. The latter is a consequence of 

' enhanced radiative damping for modes which propagate immediately below the convection zone. 

The critical value of lut c is of observational significance because modes with the smallest value of 
lot c are most observable photometrically. Maximum periods for overstable modes predicted for 
our cooler model envelopes are about a factor two longer than the observational upper limit of 
1,200 s. We assess a number of plausible resolutions for this discrepancy among which convective 
overshoot and nonlinear saturation look promising. The nonadiabatic eigenfunctions enable 
us to predict relative amplitudes and phases of photospheric variations of flux and velocity, 
quantities made accessible by recent observations. We also present asymptotic formula for 
damping rates of high order modes, a result of consequence for future investigations of nonlinear 
saturation of the amplidues of overstable modes. 
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Introduction 



While passing through a narrow temperature range around T c g ~ 12,000 K, DA (hydrogen) white 
dwarfs exhibit pulsations with periods between 2 and 20 minutes. These oscillations are identified with 
g-modes having radial orders 1 < n < 30, and spherical harmonic degrees 1 < £ < 2. Mode overstability 
occurs when driving in the thin and fast reacting surface convection zone exceeds damping in the underlying 
radiative interior (Brickhill [1991 , Gautschy et al. 1996, Goldreich & Wu 1998, hereafter Paper I). The 



flux perturbation entering the convection zone from the underlying radiative interior is in phase with the 
pressure perturbation. Because the convection is fast and efficient, the specific entropy perturbation is 
nearly depth-independent throughout the bulk of the convection zone. Consequently, the magnitude of the 
flux perturbation decreases outward. The phase relation between the pressure and flux perturbations then 
results in convective driving. 

Gravity-modes have simple structures in their upper evanescent regions. This permits an analytic 
derivation of a stability criterion based on the quasiadiabatic approximation for modes whose evanescent 
regions extend well below the base of the convection 'zone. Overstability is predicted provided ujt c > 1, 
where r c is a few times the thermal time scale at the bottom of the convection zone. Applied blindly, this 
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criterion predicts overstability in cool ZZ Cetis for modes whose periods far exceed those observed. Our 
numerical calculations enable us to drop both the quasiadiabatic approximation and the restriction that 
the evanescent region extend well below the convection zone. They also allow us to incorporate the effect of 



turbulent viscosity as discussed in Brickhill ( 1990 ) and Goldreich & Wu ( 1998 , hereafter Paper II). 



This paper is the third in a series dedicated to the overstability of gravity-modes in ZZ Cetis. It 
is largely numerical, and presents results of nonadiabatic calculations of eigenvalues and eigenfunctions. 
Nonadiabaticity is important for modes which have lot c <J 1, since they have r t h < 1 in their driving and 
damping regions. This category includes both overstable modes with wr c ~ 1, which are easily detectable by 
photometric measurements, and damped modes with ljt c <C 1, whose excitation by parametric instability 
acts to limit the amplitudes of overstable modes. It is imperative that our calculations accurately represent 
them. 

The organization of this paper is as follows. In §^|, we collect the equations which govern linear, 
nonadiabatic oscillations, establish appropriate boundary conditions, and describe the method we use 
to construct envelope models for ZZ Cetis. In §|| we assess the magnitude of nonadiabatic effects 
in the evanescent and propagating regions of g-modes. Numerical techniques employed to solve the 
eigenvalue problem are elaborated in §ji|. The impact of nonadiabaticity and turbulent viscosity on g-mode 
eigenvalues and eigenfunctions is described in §|[ We then discuss in §[| several possibilities to reconcile the 
nonadiabatically calculated spectra of overstable modes with those obesrved. A short summary follows in 
§[7]. The appendix is devoted to analysis of a toy model which elucidates the asymptotic damping rates of 
high order modes. 

Symbols used in this paper are defined in Table 1. Their usage is consistent with that in previous 
papers of this series. 



2. Equations and Boundary Conditions 
2.1. Equations in the Radiative Region 

The linearized equations describing nonadiabatic pulsations read 

'"-a = ikJ^-z z ), (2) 
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for plane-parallel geometry with the Cowling approximation, and an assumed time-dependence of exp~ luJt 
Equations (0)-(Q) express the conservation of mass, horizontal and vertical momentum, as well as energy. 
Equation rta) describes the equation of radiative diffusion. Closure of this system of equations requires 
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Table 1. DEFINITIONS 



Symbol Meaning 

R stellar radius 

g surface gravity 

z depth below photosphere 

Zb depth at bottom of convection zone 

z w depth at top of mode's cavity, z u ~ lo 2 /(gk%) 

lo radian mode frequency, u> = u) r + iui 

n radial order of mode 

I angular degree of mode 

kh horizontal wave vector, k 2 = £(£+ 1)/-R 2 

k z vertical wave vector 

p mass density, pb mass density at z = Zb 

p pressure 

T temperature 

s specific entropy in units of fcs/mj, 

F energy flux 

c s adiabatic sound speed, c 2 a — (dp/dp) s = Tip/p 

p s p s = (d\np/ds) p 

S denotes Lagrangian perturbation 

£/j horizontal component of displacement vector 

£ z vertical component of displacement vector 

Vh horizontal component of velocity vector 

/ coefficient measuring convective inefficiency 

v cv convective velocity, v cv ~ (F/p) 1 ^ 3 

t cv response time for convection, t cv ~ z/v cv 

A, B, C dimcnsionlcss constants approximately 2, 8 & 8 for ZZ Cetis 

T t h thermal constant at depth z, t cv /r t h ~ (v cv /c s ) 2 in the convection zone 

t^, thermal time at z — z^ 

Tb unconventional thermal time constant at Zb, Tb ~ r t h/5 at Zb- 

t c time constant of low pass filter for convection zone, r c = (B + C)rb 

N 2 Brunt- Vaisala, frequency 

kt, k p opacity derivatives, kt = (d\nK/dlnT) p , k p = (d\n k/OIu p)t 

c p dimensionless heat capacity, c p = {ds/d\nT) p 

V temperature gradient, V = d\nT/d\np) 

V a d adiabatic temperature gradient, V a d = (d\nT/d\np) s 



-4- 



constitutive relations for the equation of state and the opacity. Setting Ss = reduces the above equations 
to the adiabatic ones studied in Paper I. 



We choose Sp and Ss as our independent thermodynamic variables and set 

Sp 1 Sp 
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(6) 
(7) 



Then the linear perturbation equations may be written as four first-order differential equations for the four 
dependent variables Sp/p, d(Sp/p)/dlnp, SF/F, and Ss; 
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The perturbation equations are to be solved as an eigenvalue problem. Since equation ([1C|) includes a 
factor i, all four dependent variables and the eigenfrequency are complex. We write each complex variable 
in the form 

Q = (Q r +iQ l )exp- lut . (12) 



The physical perturbation is given by Re(Q) = e Uit (Q r cos LO r t + Qi smu) r t). 



2.2. Equations in the Convection Zone 

Inside the convection zone, g-mode perturbations are constrained by the rapid response of convection 



to the instantaneous pulsational state (Brickhill |1990| & |199l| , Papers I & II). Rapid momentum diffusion 
enforces 
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And fast entropy mixing ensures that 
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(13) 



(14) 



except in the thin superadiabatic layer. The entropy gradient in the superadiabatic layer is non-negligible 
and increases with increasing convective flux. But as this region is thin and has low mass, it is sufficient to 
incorporate its effect into the boundary conditions. The total, convective plus radiative, flux perturbation 
follows from equation (|l0|). whereas equation ( |ll|) determines the radiative flux perturbation in terms of Sp 
and Ss. 
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2.3. Boundary Conditions 

Solving four linear, homogeneous, first-order, ordinary differential equations to obtain eigenvalues and 
eigenfunctions, requires a total of five boundary conditions. Four of these express physical constraints 
imposed by the environment outside the domain in which the differential equations are to be integrated. 
The fifth merely sets the magnitude scale and phase for the eigenfunctions. 



Approximations described in |2.2| enable us to lower the outer boundary from the photosphere to the 
top of the radiative interior at Zb- The three boundary conditions applied there read: 

= C, (15) 
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The constant C in equation ([[5]) sets the scale and phase of the eigenfunction. Equation ( |f6| ) follows from 
the near vanishing of X in the convection zone, as expressed by equation (|l3|), together with equations (20) 
and (21) of Paper II, which account for the jump in X across the convective-radiative boundary. Equation 
( |l7| ) relates the entropy perturbation in the main part of the convection zone to the flux perturbation that 
enters from below (see Paper I). 

It is advantageous to raise the bottom boundary from the center of the star to a depth, z = Zdoop, 
where the quasiadiabatic approximation is valid and the plane-parallel approximation still applies. This 
step alters the spectrum of eigenvalues to. However, a simple procedure to be described later allows us to 
recover values for u>i appropriate to a complete stellar model. We generally take z = z^ ccp to be the level 
at which p = 10 16 dynecm~ 2 . Since r t h rs I0 10 s at Zdeep, the quasiadiabatic approximation is valid there 
for all g-modes of interest to our investigation. And in DA white dwarf models provided to us by Bradley 



(1996), the region above Zdoep extends over the outer 2 percent of the stellar radius and includes about 10 



of the stellar mass. The two boundary conditions imposed at z = z^eep are 



l- k 4) 5 * = 0, (18) 
SF ^ < M 1 (^)+M i X = 0. (19) 



F J \p 

The mechanical boundary condition given by equation (|l8| ) states that £ z = at z = Zd ccp . It follows from 
equations (j^) and (|^). Thus our model is bounded from below by a rigid wall. Equation ( |l9|) is our thermal 
boundary condition. It is the quasiadiabatic limit of the radiative diffusion equation (|ll| ) which defines the 
coefficients M\ and M^. Our lower boundary conditions are somewhat arbitrary. For example, we could 
have adopted d^h/dz — and Ss = as mechanical and thermal boundary conditions. However, because Ss 
decreases rapidly with depth, taking Ss = as the thermal boundary condition makes it more difficult for 
our numerical schemes to converge on eigenvalues. 
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Fig. 1. — Convection zone characteristics for a sequence of hydrogen envelope models covering the 
temperature range where ZZ Cetis reside. Models have surface gravity, g — 10 8 cms~ 2 , and convection 
parameter, / = 0.32. Photospheric pressure (p P h) and pressure at Zb(pb) are in units ofdynecm -2 . Time 
constants t cv , Tb, Tth, and r c are in units of s. Both t cv and rth are evaluated at Zb- Observationally detected 
g-modes have periods in the range 10 2 — 10 3 s. 



2.4. Envelope Models 

Instead of using complete white dwarf models, we work with plane-parallel, hydrogen envelopes 
computed on fine grids. These are produced by integrating downward from the photosphere^ using the 
Lawrence Livermore equation of state and opacity tables (Rogers et al. 1996 , Iglesias & Rogers |l996 ). Our 
envelopes have g = 10 s cms -2 and cover the range 11,000 K < T e ff < 13,000 K. 

We model convection by invoking the mixing length ansatz. This involves introducing a dimensionless 
parameter / such that^j 

ds / / m p F cv y ^ 

dln P {\Ps\ppY \kB T J 

The parameter / is of order unity and is related to the conventional mixing length ratio a by / ~ a~ 4 / 3 . 
A smaller / signifies more efficient convection and yields a thicker convection zone. The radiative flux is 



3 The photosphere is taken to be where p = 2g/(3n). 
4 See equation (23) of Paper I. 
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Fig. 2. — Comparisons of eigenfunctions calculated with and without the inclusion of radiative diffusion. 
The upper panel shows the fractional Lagrangian flux perturbation, 5F/F, and the lower the differential 
work, defined by W = J 5Wd\ogp, both plotted against log p. Results from quasiadiabatic calculations 
are depicted by dashed lines, whereas those from nonadiabatic calculations are given by solid lines for real 
parts and dotted lines for imaginary parts. We set the nonadiabatic 5F/F to have the same phase and 
amplitude as the quasiadiabatic one in the deep adiabatic interior. The stellar model has T e g = 12,000 K, 
which with g = 10 8 cms~ 2 and / = 0.32 yields r t h = 200 s at zj,. The mode shown here has £ = 1 
and P — 430s. It is fairly adiabatic in the radiative region since ujt^ — 2.8 at zi,. Nevertheless, radiative 
diffusion noticebly reduces convective driving and gives rise to a small region of radiative driving immediately 
below the convection zone. 



related to the entropy gradient by 

ds _ 3K P F Tad 

dlnp ~ WaT^g ad ' ( ' 

We determine ds/dlnp, F cv , and F ra( j from equations ( p0| ) and ( pi]) together with the relation F = F cv +F Ta ^. 

The grids of our model envelopes arc chosen fine enough to resolve the steep entropy profiles in the 
superadiabatic layer. Figure [l] displays some characteristics of the surface convection zones in envelope 
models produced with / = 0.32. This value of / enables us to match our model to that of Bradley's at 
T e ft = 12,420 K. Notice that the eddy turn-over time, t cv , is of order a few seconds in even the coolest 
models of interest. 
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Fig. 3. — Similar to Figure |2|, but for a mode with I = 1 and P = 1,500 s. Since ujTth — 0.8 at Zb, this 
mode is moderately nonadiabatic at the top of the radiative interior. The mode remains overstable when 
radiative diffusion is taken into account. However, the flux perturbation at the photosphere is modified in 
both amplitude and phase. 

3. Where Nonadiabaticity Is Important 

Nonadiabaticity is quantified by how much the presence of the entropy perturbation, 5s, affects Sp/p 
and 5F/F. Our analysis for the radiative interior leads to separate criteria for the evanescent region and 
the propagating cavity. 



3.1. Nonadiabatic Effects in the Evanescent Region 



The effect of Ss on Sp/p is contained in equation (|^). Noting that 5p/p varies on the scale z^ > z, we 
obtain 
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Fig. 4. — Similar to Figure g, but for a mode with I — 1 and P = 2,400 s corresponding to w r th = 0.5 at 
zj). The nonadiabatic region extends down into the propagating cavity. There is a large phase difference 
between 5F/F and 5p/p at Zb which leads to convective damping. As a consequence, this mode is damped, 
although a quasiadiabatic calculation predicts overstability. 

where z ~ p/{gp) and N 2 ~ g/z. The nonadiabatic term is of order, 

p (gk h ) 2 - cj 4 z 

p s os ~ — os. (23) 

P9 9^ Zaj 

So the nonadiabatic correction to 5p/p is of order 5s. 

To relate 5s to 5p/p, we turn to equations ([l0]) and (|Tl"|). With appropriate scalings they yield]^] 

wTth amp V ^ / WT th ^ 
f ~ (24b) 



Together, equations (|23])-( p4b| ) imply that the ratio of the nonadiabatic to adiabatic contributions to 
both 5p/p and 5F/F is of order (cjr t h) _1 , as is commonly cited in the literature. 



3 The coefficient connecting Sp/p to SF/F in equation ( pj| ) varies on scale 2, although weakly. 
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3.2. Nonadiabatic Effects in the G-Mode Cavity 



Here, all perturbation quantities vary on a vertical scale k z , where k z z > 1. Moreover, it follows from 
the local dispersion relation derived in Paper I that k z <~ (zz u ) -1 / 2 . 

Scaling the adiabatic terms in equation Q yields 
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These equations imply that the nonadiabatic correction to Sp/p is of order 5s. 
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We scale equations (10) and ( |ll[ ) to relate 5 s to 5p/p; 
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Combining equations (25c)-(27b), we determine that the ratio of nonadiabatic to adiabatic contributions 
to 5p/p and 5F/F is of order (k z z) 2 /wr t h- This is not surprising; T t h/{k z z) 2 is the timescale of thermal 
diffusion across distance kg. Nonadiabaticity is measured by the ratio of the mode period to this timescale. 



4. Solution Of Eigenvalue Problem 

We follow a two step procedure in solving the linear pulsation equations (|8|)-(pd|) for the four dependent 
variables subject to the five boundary conditions given by equations (]l5|)-(]l9|). The initial step consists 
of guessing a value for the complex eigenfrequency and then applying a relaxation method (cf. Press et 
al. 1992| ) to solve the differential equations subject to four out of the five boundary conditions.^ The step is 
complete when the dependent variables satisfy both the pulsation equations and the boundary conditions 
to 10~ 7 of a scaling factor provided by the corresponding adiabatic eigenfunction at individual points. 
Working in double precision, this is routinely achieved. The second step determines the eigenvalue by 
requiring the remaining boundary condition to be satisfied. Both minimization and root finding techniques 
work well. Normally we reserve equation ( |l8| ) for our fifth boundary condition, but identical results are 
obtained when others are used instead. 



3 The trivial boundary condition given by equation ( |l^ ) is always included as one of this foursome. 
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4.1. Eigenvalues for Stellar G-Modes 

Nonadiabatic eigenvalues of our plane-parallel hydrogen envelopes are denoted by primes, a/ = Lu'j. + iuj^. 
To relate u>[ to oji of a complete DA white dwarf model, we proceed as follows. We compute adiabatic modes 
for a corresponding complete model and indentify those whose frequencies bracket oj' r ^\ We determine by 
interpolation the non- integer number n of nodes in Sp/p to assign to the adiabatic mode of the complete 
model whose frequency uj — ui' r . Suppose that n! of these nodes lie above Zdccp- Then 

n! 

Ui « — w{. (28) 



Validation of equation feq) is provided by equation ( A7) in the appendix. Its physical justification is that 



whereas energy is dissipated near the surface, it is stored throughout the mode cavity. 

Our procedure for converting oj[ to u>i requires nonadiabatic effects to be small well above Zdccp, as 
is the case for all modes of interest to our investigation. However, it only applies to modes which have at 
least one node above Zdccp- Thus it cannot handle the highest frequency g-modes. Values of u>i for these are 
computed from the work integral. 

Very low order modes (n < 6) are evanescent above Zdeep- We therefore cannot enforce the boundary 
condition equation on them to solve for the eigenvalues. We rely on work integral for stability analysis 
of these modes. 



4.2. The Work Integral 



Calculating the work integral provides an approximate method for evaluating the driving or damping 
rate of a mode. The work integral computes 7 = 2u>i. This method is well defined and accurate when the 
quasiadiabatic approximation applies. Moreover, it reveals regions of driving and damping. Since overstable 
g-modes of ZZ Cetis have high quality factors, their linear pulsations are nearly periodic. The work integral 



may be calculated as follows (Unno et al. |1989[ ) 
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5. Numerical Results 

5.1. Nonadiabatic Effects On Eigenfunctions 

The effects of nonadiabaticity are illustrated by comparing nonadiabatic and adiabatic eigenfunctions 
for three 1=1 g-modes of the same stellar model. The model is characterized by T c g = 12,000 K, 
g = 10 8 cms~ 2 , and / = 0.32, which together imply r t h ~ 200 s at z = z&. The modes have periods of 430 s, 
1,500 s, and 2,400 s. Quasiadiabatic calculations predict overstability for each of these modes. 



7 The adiabatic modes have the same I as the mode of the plane-parallel envelope. 
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Fig. 5. — Amplitudes and phases of 8F/F and Vh as functions of logp. The displayed mode has I = 1 and 
P = 800 s corresponding to wrth = 0.8 at z = z;, in a stellar model with T e g = 12,160 K, g = 10 s cm s~ 2 , and 
/ = 0.32. The amplitude, \X\, and phase, <fi, of a perturbation X are defined as X = \X\ cos(ujt — </>). In the 
adiabatic interior, the flux perturbation leads the velocity perturbation by 90°. 



Radiative diffusion acts to soften sharp temperature gradient perturbations. Its importance increases 
with mode period. These characteristics are illustrated in Figures [|-|}f] The 430 s mode, with uiT t h — 2.8 at 
z = zi,, is quite adiabatic in the radiative interior. Radiative diffusion is more pronounced in the evanescent 
region of the 1,500 s mode, which has UTth = 0.8 at z = Zb- The nonadiabatic eigenfunction of the 2,400 s 
mode, for which Lur t h = 0.5 at z = Zb, deviates significantly from its adiabatic counterpart in both the 
evanescent and propagating regions. The phase shift between SF/F and Sp/p at 2% is so large that the 
convection zone contributes to mode damping. This mode is stable. 



The recent detection of velocity signals associated with g-modes by van Kerkwijk et al. fll998| ) is an 



important advance in the study of ZZ Cetis. What is being detected are horizontal velocities near the 
stellar limb. Phase differences and amplitude ratios of flux and velocity variations due to individual modes 
offer unique clues about both the star and the modes. Nonadiabatic calculations, such as that depicted in 
Figure |[ provide a theoretical basis for interpreting these observable quantities. Deep inside the radiative 
interior, the maximum flux perturbation preceeds the maximum horizontal velocity by a quarter of a 
period, and the flux to velocity amplitude ratio is determined by the quasiadiabatic approximation. These 



8 We normalize Sp/p at z = Zd eep by setting it equal to its adiabatic counterpart. Thus the imaginary component of 6F/F 
is due to nonadiabaticity. 
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relations are modified in the upper envelope by the convective transport of heat and momentum and by 
radiative diffusion. Convection delays and diminishes the flux variation (see Paper I). It also flattens the 
horizontal velocity profile above zj, and forces it to jump at Zf, (see Paper II). Radiative diffusion smears 
flux perturbations below the convection zone. In the example illustrated in Figure ^, a photospheric flux 
variation of 20mma[] corresponds to a horizontal velocity of Vh — 2.3 km s -1 . Nonadiabaticity, mostly 
in the convection zone, reduces the phase lag between velocity maximum and light maximum from the 
quasiadiabatic value of 90° to 55° at the photosphere. Employing the conventions in van Kerkwijk et 
al. ( |1998| ), this mode would exhibit R v = \v h \/(uj\5F/F\) = 16, and A$ v = $(8F/F) - $(«/,) = -55°. This 
is consistent with the 818 s mode which they observed to have Rv = 11 ± 4 and A$„ = —44° ± 19° .[^] 



5.2. Nonadiabatic Effects On Driving and Damping Rates 

Where applicable, direct calculations of u>i yield results consistent with those based on the work 
integral. This provides a measure of confidence in both. Modest discrepancies are found for some marginally 
overstable modes. The eigenvalue calculations are more reliable, as work integrals suffer from inaccuracies 
due to cancellation between comparable magnitudes of driving and damping. Values of LUi for individual 
g- modes evolve with decreasing T c r. Effects of nonadiabaticity on driving and damping rates of g- modes 
are presented in Figures |^-^. General trends are described below. 

Upper lids of cavities of short period modes (e.g., the 430 s mode) lie far beneath the bottom of the 
convection zone in even in the coolest variables. Their driving/damping rates are largely immune to both 
nonadiabatic effects and the depth of the convection zone, and hardly vary across the instability strip. The 
quasiadiabatic estimate, u)\ ~ pertains to these modes (see Paper I). 

Longer period modes (e.g., the 1,000 s mode) become overstable at lower T e g. Values of their \ui\ 
exhibit a steady decline with decreasing T e g. This is a consequence of the increase in mode inertia with 
decreasing T e g . As the convection zone deepens, it depresses the upper lid of a mode's cavity. This decreases 
the relative size of the perturbation amplitude near the surface, where nonadiabaticity is greatest, with 
respect to that in the interior, where most of the mode energy is stored. 

Modes with even longer periods (e.g., the 2,000s mode) are weakly overstable in cooler and narrower 
temperature ranges. Nonadiabatic effects tend to stabilize these modes. In the limit of strong dissipation, 
u>i/uj r w — \nlZ~ 1 / (2irn), where 1Z is the reflection coefficient at the cavity lid. We derive this relation 
using a toy model in the appendix. Magnitudes of nonadiabatic damping rates are an important input 
to calculations of parametric instability, an amplitude limiting mechanism for overstable modes that is 
explored in the next paper of this series. 

6. Maximum Periods of Overstable Modes 



9 Milli-magnitude-amplitude is the commonly adopted unit for measuring flux perturbations. An 1 mma variation means 
8F/F = 10" 3 . 

10 Effects of limb-darkening, disc averaging, and bolometric correction introduce a factor close to unity. 
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Fig. 6. — Driving and damping rates for 1=1 modes plotted against mode periods. Top to bottom, the 
panels pertain to stellar models computed with g = 10 8 cms and / = 0.32 for T c s = 12,600 K, 12,300 K, 
and 11,750 K. Values of \oJi\ obtained from the nonadiabatic code are denoted by solid triangles. Those 
calculated from work integrals using nonadiabatic eigenfunctions are shown by open squares. For very short 
periods, P < 400 s, our shallow envelopes force us to rely entirely on the work integral. The dashed vertical 
line marks the boundary between shorter period overstable modes and longer period damped ones. The 
solid line displays the analytic estimate, cj, = — uj r hxli," /(27m), for the damping rate of modes which suffer 
strong dissipation. We set 1Z — 1/e here, although in reality it decreases with increasing mode period and 
effective temperature. 



6.1. Limitations of the Overstability Criterion lot c > 1 
The derivation of the overstability criterion lut c > 1 depends on the quasiadiabatic approximation and 



the condition that z u 3> z\, (Brickhill 1991, Paper I). Since r c w 4r t h at Zb, the quasiadiabatic approximation 
is suspect for modes with lut c w 1. The upper and lower panels of Figure ^ display the boundary between 
overstable and damped modes as determined by quasiadiabatic and nonadiabatic calculations. Both show 
that the boundary value of lot c is close to unity for hot ZZ Cetis but rises as stars cool. Inaccuracy of the 
quasiadiabatic approximation does not dramatically modify the overstability criterion, but the violation of 
the condition z u 3> z& does. The latter occurs because Zb increases with decreasing T c g and z^ decreases 
with decreasing oj. Modes with z u < Zb propagate just below Zb- Their short WKB wavelengths enhance 
radiative damping above the estimate given by equation (54) in Paper I on which the derivation of the 
overstability criterion lut c > 1 rests. 
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Fig. 7. — Panoply of overstable I = 1 modes according to quasiadiabatic (upper panel) and fully nonadiabatic 
(lower panel) calculations in white dwarf models produced with g — 10 s cms -1 and / = 0.32 for a range of 
T e ff which covers the ZZ Ceti instability strip. Each overstable mode is marked by a solid triangle whose size 
is proportional to the diminution factor of the flux perturbation in the convection zone, [1 + (wtc) 2 ]" 1 / 2 (see 
Paper I). Dashed lines denote the loci of constant values of lot c . The dotted line corresponds to ujt^ = 1 
with Tth evaluated at z = Zb- 



6.2. Comparison with Observations 



Both our quasiadiabatic and nonadiabatic calculations consistently predict maximum periods for 
overstable modes that are about a factor two longer than those observed (cf. Fig. ||). We have not been 
able to resolve this discrepancy. Several possible explanations which we have considered are described in 
what follows. 



6.2.1. Effects of Turbulent Damping 



Both Brickhill (199C) and Paper II stress that linear damping due to turbulent viscosity, although 
negligible inside the convection zone for all modes, may be significant in the region of convective overshoot 
for long period modes. We incorporate this effect into our analysis, while recognizing that it cannot be 
quantified precisely. 
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Fig. 8. — Observational and theoretical instability strips for DA white dwarfs. The upper panel displays 



each of the known ZZ Cetis according to the dominant period in its pulsation spectrum (Bradley 1996), 
and its inferred effective temperature (Bergeron et al. 1995). The lower panel repeats the material shown 
in the lower panel of Figure ^. The dotted and dashed lines illustrate the maximum period for overstable 
modes under the stabilizing effect of turbulent damping in the region of convective overshoot for values of the 
parameter X/zb equal to 0.1 and 1.0, respectively. The dotted-dashed line shows the envelope of overstable 
modes for models which incorporate one pressure scale height of convective overshoot. When comparing 
observational and theoretical instability strips, one should bear in mind that both depend upon the assumed 
mixing-length. 



The rate of turbulent damping in the overshoot region is expressed as (eq. [36] in Paper II), 



7v 



-wR 2 p b UJ 3 X\A^ h \ 



(30) 



where |A^| is the normalized jump in the horizontal displacement across the convective overshoot region 
within which the turbulent viscosity decays on length scale A. To cover our ignorance, we present results 
for X/zi, of 0.1 and 1.0. Figure || demonstrates that with the larger value, the longest period overstable 
mode at a given T e g is compatible with the observed upper limit of P — 1,200 s. However, it seems unlikely 
that X/zi, w 1. The overshoot region probably extends less than a pressure scale height below z^, while 
Hp w OAzb, and A must be several times smaller than H p . 
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Fig. 9. — Effects of convective overshoot on g-mode structure and stability. The stellar model lies at 
T c ff = 12,250 K with one pressure scale height of overshoot below the convection zone. We use the same 
mode as in Figure 0, and the stellar model there has a convection zone as deep as the bottom of the overshoot 
region in this model. The small inset in the top panel shows the Brunt- Vaisala frequencies for the former 
in dashed line (labelled with 'b') and the latter in solid line (labelled with 'a'). The top and bottom panels 
display similar materials as in Figure except that the nonadiabatic work integral is multiplied by a factor 
of 20 to aid display. Both quasiadiabatic and nonadiabatic calculations predict stability for this mode. 



6.2.2. Effects of Convective Overshoot 

We modify our envelope models to simulate convective overshoot by adding an isentropic layer below 
the base of the convection zone. Within this layer the radiative flux exceeds the total flux and turbulence 
transports energy downward. The top of the radiative interior immediately below the overshoot layer 
differs from that below a convective layer of similar depth. The temperature gradient is shallower and 
the Brunt- Vaisala frequency has a finite positive value. These differences enhance the effects of radiative 
diffusion and modify the relations between horizontal velocity and light perturbations at the photosphere. 

Figure ^ provides such an example. Our nonadiabatic calculations indicate that one pressure scale 
height of overshoot suffices to stabilize all modes with P > 1,400 s (see Fig. |J). Moreover, by effectively 
deepening the convection zone, it shifts the instability strip to higher T e g. 
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6.2.3. Sensitivity to Surface Gravity 

Our nonadiabatic calculations indicate that, at fixed T c g, the longest period for an overstable mode 
scales approximately as g^ 1 ^ 3 . Thus maximum periods of overstable modes around 1,200 s would require 
g w 5.0 x 10 s cms -2 , which is well outside observational constraints.^] 



6.2.4- Other Observational and Theoretical Considerations 



Hansen et al. (1985) show that damping due to upward propagation of a gravity wave above the 
photosphere is unimportant for modes having frequencies exhibited ZZ Cetis. We concur with this 
conclusion. 

Detection of lower frequency modes requires longer observational data streams. Noise due to variations 



of atmospheric transparency increases at lower frequencies (see Fig. 4 of Winget 1991 ). These may result 
in an observational bias against the detection of low frequency modes. 

Evidence that photometric amplitudes of modes decline for P J> 1,000 s, suggests that the observational 
cutoff at P pa 1,200 s) is genuine (see Fig. 5 in Clemens 1995| ). Taken literally, it also hints that the cutoff 
is the result of a nonlinear mechanism which saturates pulsation amplitudes. 



7. Summary 



There is little doubt that convective driving, as originally proposed by Brickhill (1991), is responsible 
for the linear overstability of g-modes in DA white dwarfs. It is physically self-consistent, and convincingly 
rationalizes observational facts. It accounts for the general location of the instability strip, although 
a precise specification depends upon the modeling of convection (e.g., the mixing-length parameter). 
Convective driving also explains why longer period modes become overstable as a star cools (see Fig. |^). 

Although we agree with Brickhill that ojt c > 1 is a necessary condition for mode overstability, we find 
that it is not a sufficient condition for modes whose periods exceed 1,000 s. This stems from enhanced 
radiative damping of modes whose upper cavity lids approach as is apparent from both quasiadiabatic 
and nonadiabatic calculations (see Fig. [?]). Our nonadiabatic calculations yield a maximum period of about 
2,300 s for overstable modes. This clashes with the maximum period of 1,200 s for observationally detected 
modes. 



We also agree with Brickhill's deduction (1990) that turbulent convection forces the horizontal velocity 
to be nearly independent of depth within the convective envelope. Consequently, mode damping due to 
turbulent dissipation within the convection zone is reduced to a negligible level. However, suppression of 
the horizontal shear in the convective envelope results in a shear layer at the top of the radiative interior. 
When convective overshoot is accounted for, this provides linear turbulent damping. Figure || suggests that 
turbulent dissipation might depress the maximum period of observable modes to some degree. Convective 
overshoot also alters the thermal structure of the upper radiative layer. A rough treatment of this effect 
with one scale height of overshoot predicts a maximum period of 1,400 s for overstable modes. Nonlinear 



11 Applying a theoretical mass-radius relation (Hamada & Salpetcr 1961) for white dwarfs, this corresponds to a stellar mass 
about 1.2M . 
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interactions which limit the amplitudes of overstable modes may also play a part in determining the 
maximum period. 

In conclusion, the blue edge of the theoretical instability strip seems to be set by the condition 
that uot c w 1 for the lowest order I = 1 mode. But the location of the red edge is more nebulous, and 
may result from a combination of decreased photometric visibility, convective overshoot, and nonlinear 
effects. Additional detections of velocity variations associated with g-modes could provide important clues. 
Convective driving makes the testable prediction that velocity variations become relatively more observable 
than photometric variations toward the red edge of the instability strip. 
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A. A Toy Model For Nonadiabatic Modes 

We describe a simple toy model for nonadiabatic modes. It is particularly useful for interpreting 
damping rates in the limit of strong dissipation. 

Consider waves which satisfy the one-dimensional, homogeneous, acoustic wave equation 

in the interval < z < L. Here £ is the Lagrangian displacement and u is the constant propagation speed. 
The dispersion relation connecting frequency, w, and wave vector, k, reads to 2 = k 2 u 2 . We take the lower 
boundary to be a rigid, perfectly reflecting wall, so 

£ = at z = L. (A2) 

Dissipation is introduced by means of a partially reflective upper boundary, where 1Z denotes the reflection 
coefficient of the incident wave. This is expressed through the boundary condition 

d 3\ . „ ( d d 



Note that in the limit 1Z —> 1 , the upper boundary becomes a perfectly reflecting wall. 
Eigen-solutions of equation ( |Al[ ) are composed of oppositely directed waves; 

£ = j± e -iut-ikz _|_ g e -iut+ikz (A4) 



Application of the boundary conditions given by equations (A2) and (A3) yields 

B = -HA, (A5) 

and 

kL = nix - ^lnft -1 , (A6) 
where n is the number of half wavelengths between the walls. Then the dispersion relation implies 

lo 1 = -±u r = lnT^ -1 , (A7) 

K r Zirn 

where n is an integer. 



Equation (A7) is the key result of our toy model. It demonstrates that u>i grows logarithmically with 
7^ _1 , and provides an order-of-magnitude estimate for the damping rate of strongly nonadiabatic stellar 
modes. Radiative diffusion substantially suppresses the effective reflection coefficients of modes having 
ut^j < 1. In hot DA stars, this applies to modes with periods in excess of 1,000 s (cf. Fig. ||). Finally, 
equation (|A7| ) justifies the procedure we apply in §4T to translate the damping rate for a mode of an 
envelope model, to u>i, the damping rate for a mode of a complete stellar model. 



